Thermodynamics and phase behavior of the lamellar Zwanzig model 
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Binary mixtures of lamellar colloids represented by hard platelets are studied within a general- 
ization of the Zwanzig model for rods, whereby the square cuboids can take only three orientations 
along the x, y ox z axes. The free energy is calculated within Rosenfeld's "Fundamental Measure 
Theory" (FMT) adapted to the present model. In the one-component limit, the model exhibits the 
expected isotropic to nematic phase transition, which narrows as the aspect ratio C, = L/D {D is the 
width and L the thickness of the platelets) increases. In the binary case the competition between 
nematic ordering and depletion-induced segregation leads to rich phase behaviour. 

PACS numbers: 61.20.-p, 61.30.Gd, 82.70.Dd 



I. INTRODUCTION 

, Suspensions of lamellar colloids, like clay dispersions, 
'have received renewed experimental and theoretical atten- 
|tion lately, because of the expected rich phase behaviour, 
■and the obvious geophysical jand technological implications, 
in particular for oil drillingtl With increasing concentra- 
tion moat natural and synthetic clay dispersions undergo 
gelationo which impedes the isotropic to nematiG-|(I-N) tran- 
sition expected ffffln theoretical considerationa3 and com- 
puter simulationsQEI on entropic grounds. It came hence as a 
relief when a novel system of hard lamellar colloids, namely 
■a dispersion of gibbsite platelets sterically stabilized by a 
^grafted polymer layer, in toluene, was shown to exhibit the 
expected I-N transitionQ. The difference in behaviour be- 
tween the gel-forming clays and the novel gibbsite system 
jnay lie in the long-range Coulomb interactions present in 
■the aqueous clay dispersionglL Subsequently it was-shown 
|experimentallyl3, theoreticallyu, and by simulationll3, that 
polydispersity in the size of the platelets strongly affects the 
phase behaviour. In particular binary mixtures of thin and 
.thick platelets, kad to an unexpected isotropic-nematic den- 
'sity inversionH'Q. Moreover, even for monodisperse systems, 
the phase diagram was shown to be quantitatively, and even 
qualitatively, sensitive to the aspect ratio (i.e., the thickness 
pver diameter of platelets modeled as " cut spheres" )eI. 
' In this paper we address the problem of the phase be- 
Jiaviour of monodisperse and bidisperse systems of hard 
platelets using a highly simplified model introduced some 
'time ago by Zwarpig to study the I-N transition in sys- 
tems of hard rodal^. In this model platelets or rods are 
represented by square parallelepipeds (or cuboids) of thick- 
ness (length) L and width D, which can take only three 
orientations, along the x, y or z directions, rather than a 
continuous range of orientations in space (cf. Fig.0). The 
aspect ratio ( = L/D > 1 for rods, while C < 1 for platelets. 
The case C 1 corresponds to a system of parallel cubes, 
for which analyti|C-expressions of the first seven virial coeffi- 
cients are knownE3. Zwanzig's model may be considered as 
a highly coarse-grained version of the Onsager long-rod sys- 
tem with continuously varying orientation. The advantage 



FIG. 1: In the model under consideration platelets are repre- 
sented by square parallelepipeds of thickness L and width D, 
which can take only three orientations, along the x, y or z direc- 
tions. The number density of the centre of mass of the platelet 
displayed in the figure is denoted by py. 



is that virial coefficients beyond the second coefficient B2 
can be calculated in the limit C ^ 1 thus providing a quan- 
titative -test of Onsager's classical theory, which includes 
only B2U. In this paper we adapt Zwanzig's model to the 
platelet case (C < 1) to investigate the equation of state, 
the TN transition, as well as binary mixtures of platelets 
of different size in the search of a possible, depletion-driven 
demixing transitiontj. 



II. MODEL AND VIRIAL EXPANSION 

Following up on Zwanzig's original idea for rods, we con- 
sider a system of hard square cuboids, of dimension 
D X D X L, as shown in Fig. |l|, with i < D (C < 1). 
The normals to these platelets are restricted to point along 
one of the Cartesian axes, and platelets of a given orienta- 
tion are considered as belonging to one of three "species" ; 
the numbers of platelets belonging to each species are de- 
noted by Nx , Ny and , with the total number of platelets 
equal to N = + Ny N^. N is fixed, but the Na 
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fluctuate due to thermal reorientations of the platelets. 
The mean numbers of platelets per unit volume will be 
denoted by —< Np > /V. In an isotropic phase, 
Px ^ Py ^ Pz ^ p/3, where p = N/V. 

The virial expansion of the excess free energy per unit 
volume, in powers of the partial densities pp, is: 



VksT 



oo ^ 

^—jBi{(3i,...,Pi)pp^ ■■■ Pl3,, 

i=2 



(1) 



where the indices /32, ■■■ sum over x, y, and z. _B2(/3i, /?2) 
is the second virial coefficient between platelets of species 
/3i and 132'- 



(2) 



where for cuboids, the Mayer functions reduce to products 
of Heaviside step functions: 



and the resulting equation of state reads: 
P 



p + B2{x,x) (pI+pI + pI) 

+ B2{X, y) {pxPy + PxPz + PyPz) 

+B3{x, X, x) {pI + pI + pI) + -83 (x, y, z)p^pyp^ 

+ B3{x,X,y) [pliPy + Pz) + pliPx + Pz) + pUPx + Py)] 

+Oip'). (7) 

It is instructive to consider the isotropic (low density) phase 
first, where Px — Py — Pz = p/3. Introducing the total (or 
orientation-averaged) virial coefficients: 



B2 = ^[B2{x,x)+2B2{x,y)] , 



(8) 



S3 = -[B3{x,x,x)+QB3{x,x,y) + 2B:i{x,y,z)] , (9) 

the equation of state of the isotropic phase reduces to the 
usual form: 
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(3) 



In Eq. (^), Sa^p is the spatial extent of plates of species 
(orientation) (3 in the direction a = x,y, z: 



(4) 



while \ra\ is the absolute value of the a-component of the 
vector Fq,. For obvious symmetry reasons, there are only 
two independent i?2(/3i, /32), namely B2{x,x) = B2{y,y) = 
B2{z,z) (parallel platelets) and B2{x,y) = B2{x,z) — 
B2{y,z) (orthogonal platelets). The corresponding integrals 
(H) are easily calculated, to yield: 



B2ix,x)^ALD^; B2{x,y) ^ D{L + Df 



(5) 



Similarly, there are 3 independent third virial coefficients, 
namely B^^x, x, x), B^^x^x^y) and B^^x^y, z). These are 
calculated in Appendix A, together with independent fourth 
virial coefficients (the latter only in the limits L — > 0, D — > 
0, and D — L). Gathering results, the reduced excess free 
energy density may be cast in the form: 
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B2{x,x) {pI + pI + pD 

+B2ix, y) [pxPy + PxPz + PyPz) 

+ ^Bz{x, X, x) {pI + pI + pI) + ^Bsix, y, z)p^pyp^ 

+ ^-B3(x,X,?/) [pliPy + Pz) + pliPx + Pz) + pUPx + Py)] 



-0{P'): 



(6) 
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p + B2P^ + Bsp^ + 0(p4) , 



(10) 



The two important dimensionless parameters in the problem 
are the size ratio ^ = L/D, and the reduced density (or 
packing fraction) p* = pD^L. Substituting the results from 
Appendix A, the dimensionless equation of state P/ {ksTp) 
may be expressed in terms of these two variables as: 



P 



ksTp 
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l + ^(i+4 + C)/ 
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141 + 52C ) +0(/^) 



(11) 



It is worthwhile to examine three limiting cases, including 
contributions up to the fourth virial coefficient. 

a) Long rods, corresponding to the limit ^ — > 00; in that 
case the relevant variable is the Onsager parameter, or ef- 
fective packing fraction 77 = p*C, = pDL^, in terms of which: 



P( 



rods) 



ksTp 



= 1 



2 16 
3^ ^ 243^ 



(12) 



i.e., the third virial contribution ^ if is identically zero, in 
agreement with jOnsager's conjecture for rods-.with continu- 
ous orientationsEl, and with Zwanzig's resultlld. 

b) Thin platelets, corresponding to the limit ^ — > 0; in 
that case the relevant variable is the effective packing frac- 
tion rj = p* /C ~ pD^, 

p(platclcts) 9 9 09 

In this case the third virial coefficient B^ is non-zero, nor is 
-B4, in agreement with numerical findings for continuously 
rotating disc-shaped platelet^. At this stage there is no 
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FIG. 2: Low density equation of state of a monodisperse, 
isotropic fluid consisting of rectangular platelets of surface size 
D X D and thickness L for various aspect ratios: C, = LjD = 
0,0.3,0.6, 1 (solid lines from bottom to top). The symbols rep- 
resent cpijiputer simulation data for parallel hard cubes liejrcles, 
73 — L)Efl and thin circular platelets (squares, L 0)u'u with 
continuous allowed orientations. 



indication of convergence of the virial series at high pack- 
ing fraction, but some plausibility arguments will be put 
forward in the following section to justify truncation of the 
virial series ( p^ ) after the 0{rj^) (third virial) term for in- 
finitely thin platelets. 

c) Parallel cubes, corresponding to C = 1; the relevant 
dimensionless variable is the packing fraction rj = p*, and 
the first few terms in the virial series readtHI: 



framework of Rosenfeld's FjMT0, by generalizing Cuesta's 
work for hard parallel cubesEII to the case of Zwanzig's model 
for hard platelets. The second and third virial coefficients 
calculated in the preceding section provide key input into 
the theory. Since the free energy of bidisperse platelets will 
be required to construct complete phase diagrams in sec- 
tion V, the FMT is formulated here for a binary mixture 
of platelets of size Di x Di x Li {i = 1,2). It is conve- 
nient to consider the more general case of inhomogeneous 
mixtures, characterized by the local densities p^*^ (r) of the 
centres of mass of platelets of species i, with orientations 
along (3 — x,y, z. 

Under the influence of external potentials vjf'^r), the 
equilibrium density profile of the mixture minimize the 
grand potential functional: 



fi[p«(r)] = / dvpf{v)[kBT(\n[klpf{v)]~l 

1=1 fj '' 



(15) 



where the pi are the chemical potentials of the two species 
(irrespective of their orientation), while the are the ther- 
mal de Broglie wavelengths. The basic assumption of FMT 
is to postulate the following form of the excess free energy 
functional: 



F,^[pf{v)]^kBT J dr<P[ni{r)]. 



(16) 



where the reduced excess free energy density $ is a function 
of a set of weighted densities 



ni 



W = EE/ dnpf{r,)u;l${r - n) . (17) 
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bcs) 
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1 + 477 + 9?7^ + —rf + 0{rf) . (14) The weight functions Lul^Ur) are determined by expressing 
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the Fourier transforms of the Mayer functions f^pfp.^ as sums 
of products of functions associated with single platelets. In 
the case of mixtures, equation (||) is generalized as: 

4:l('-i>r,) = 

Y{^{\{s'X+s'^!p)-K,,~r^a\) , (18) 



In Figure 2 the low density equation of state is shown 
for several aspect ratios < C < l,r-^nd compared to avail- 
able simulation data for hard cubeslia and for infinitely thin 
discs of the same ar|e^ {ttR^ = D^), but, with continuously 
varying orientationsQLl. As expected, the pressure increases 
rapidly with ^ for any reduced density pD^. The results 
for continuously rotating discs and for squares with discrete 
orientations are surprisingly close. 



III. FUNDAMENTAL MEASURE THEORY FOR 
HARD PLATELETS 

In view of the uncertain convergence of the virial series 
of the equation of state of hard platelets at high densities 
{v ^ l)j it is important to seek an approximate resumma- 
tion of |the series, similar to scaled particle theory for hard 
spheresEj. This is most conveniently achieved within the 



where S'^'-|g = Di + {Li — Di)5ajj- The Fourier transform of 
the Mayer function ( p^ ) may be decomposed according to: 

4:i(q) = /^---/£i(r) 

—00 

+^%, (q) • (q) + (q) • (q) ^ m 
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where the scalar and vectorial weights (0 < Z < 3) are de- 
fined by: 



(20) 



(21) 



/3fe)a,,0(g^) / 



and: 



(22) 



(23) 



(24) 



(25) 



Substitution of the Fourier transforms of the weight func- 
tions (p0|)-(|2^) into equation ( p7| ) yields two scalar and two 
vectorial weighted densities niij). Letting n/(r) = n;(r) • u, 
I = 1, 2, with u = (1, 1, 1), one can construct the function <i> 
from the following linear combination of scalar terms, each 
of dimension (volume)^^: 

$[n;(r)] = 

cono(r) + cini(r)n2(r) + C2ni(r) • n2(r) + C3n^(r) 
+C4n2n2{r) ■ n2(r) + C5n2(r) • n2(r) • n2(r) , (26) 

where n-n-n — n^+Uy+nl. The 6 coefficients Q (0 < i < 5) 
are analytical functions of the only dimensionless weighted 
density, namely 71,3. These functions are determined, within 
integration constants, by the scaled particle conditioijlid, 
which results in the following differential equation for $113: 



(27) 
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Substitution of ( |26| ) into (27) yields, upon integration: 

di d2 



Co 

C3 



ln(l - ns) ,ci 



, C4 



1-713 
(1-^3)2 



,C2 



,C5 



1 - 713 ' 

(l-na)- 



.(28) 



The remaining constants (1 < i < 5) are determined 
by identifying the low density expansion of the excess free 



energy, defined by Eqs. (|16|), ( p6D and ( |28| ) for the homoge- 
neous fluid (i.e., Pp\v) = p^^^) with the third order virial 
expansion (|l]) of the same free energy; this leads to 



di = , d2 = 1 , ^3 = ^ , ^4 = - ^ , ^5 = ■ 
D 2 3 

The final expression for the function $ then reads: 

ni • n2 n2,xn2,yn2,z 



-7io ln(l - 713) + 



1 - "3 (1 - n^Y 



(29) 



(30) 



where 112,0, is the projection of the vector n2 in the a direc- 
tion. The FMT expression for the free energy functional is 
now completely specified. In the limit of parallel hard cubes 
[Li = Di) th&,excess free energy functional reduces to the 
result of Ref.ll3. 

In the limit of a homogeneous bulk fluid the equilibrium 
density profiles are constant, and the weighted densities take 
the explicit expressions: 



(31) 



"0 = Y.{p'i^+Py 

i=i 


+p^') 
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D^P^ - 


f Api' 


2 

i=i 




- LiPy ' - 


f Api' 








+ L,pi" 






-upi^- 


f Api'' 


2 

i=l 


Upi^ + 








V Uf^ + 


^iPy " 





(32) 



p^^+p« 



(33) 



(34) 



The weight functions may be interpreted in terms of the fol- 
lowing geometric characteristics of the particles: their mean 
diameter {Li+2Di)/3, their surface 2Di{Di+2Li), and their 
volume LiD^. 

In the limit of an isotropic [p^ — Py — Pz — p/?>) one- 
component system of platelets, the equation-of-state de- 
rived from the free energy density ( ^ ) takes the following 
form (with 77 = pD^ , the effective packing fraction, and 
Q = L/D): 



P* 



= V 



(2 + new 

^ 1 - 7?C ^ 3(1 - rjcy 
2(1-^20^773 

27(1-77C)3 ■ 



(35) 



Interestingly, in the C ^ limit (infinitely thin platelets) the 
FMT equation-of-state reduces to the virial series (|l3|) , trun- 
cated after third order. This is reminiscent of the equation 



of state derived for circular platelets (discs), with continu- 
ously varying orientation, from PRISM (polymer reference 
interaction site model), which also reduces to a three term 
series (although the coefficients of rj^ and rf diSor slightly 
from the exact virial coefficient for that case)ta. These 
two observations prompt the conjecture that a three term 
virial series may be quantitatively accurate for infinitely thin 
platelets (with discrete or continuous orientations) beyond 
the low density regime. 

For finite thickness, however (^ > 0), the three-term virial 
series and the FMT equation-of-statc (^5|) deviate increas- 
ingly as C and 77 increase. Thus, for 77 = 1, Pj*, as calculated 
from Eq. ( pi] ) and i derived from Eq. (^5|), take the values 
Pi* = 2.246 and P| = 2.315 for ( = 0.1, while Pj* = 2.88 and 
Pj"^ = 3.209 for C = 0.2; for C = 0.5, the three-term virial 
series underestimates the pressure by more than a factor of 
2 (P* = 5.676 and P2* = 12.74) for ?/ = 1. 



IV. NEMATIC ORDERING OF PLATELETS 



As a direct application of the FMT free energy (30) de- 
rived in the preceding section, we now consider the possibil- 
ity of an isotopic to nematic (TN) transition in a monodis- 
perse fluid of platelets, as a function of the the aspect ratio 
C, = L/ D. The broken symmetry in the nematic phase may 
be characterized by the order parameter: 



(Pa; +Py)/2 



(36) 



which varies between (isotropic phase) and 1 (when all 
platelets are oriented along the 2;-axis). In the nematic 
phase Px = Py, and the independent variables are p and 
s. For given values of p, the one-component version of the 
grand potential 



13) (with V/5(r) = 0, and pp{r) = pp and 
given by (|16) and (^0|)) is minimized with respect to 
the order parameter s. If the FMT excess free energy ( pO[ ) 
is replaced by the virial expansion (|^) , truncated after third 
order (which is a good approximation for thin platelets, 
^ <C 1), the following Euler-Langrange equation determines 
s as a function of p: 

l + 2s 



In- 



1 - s 



2D{L-Dfps 



(D + 8L){L - DYp^s 



{L - DYp's' . (37) 



For low densities p, only the isotropic phase (s = 0) is sta- 
ble. At a critical density Pc^ a second, non-zero solution 
of Eq. (|3^) emerges, which eventually will lead to a stable 
nematic phase, when the free energy associated with s 7^ 
drops below that of the isotropic phase. The non-zero root 
s requires a numerical solution of Eq. (^7|). A simple esti- 
mate of the critical density required to obtain a non-zero 
value of s is obtained from the small s expansion of the l.h.s 
of Eq. (p7[). If only terms to linear order in p and s are 
retained, the bifurcation density is found to be: 

3 



2D{L-D) 



(38) 




FIG. 3: Equation of state of a monodisperse fluid consisting of 
rectangular platelets of surface size D x D and thickness L for 
various aspect ratios: — L/D = 0.01, 0.1, 0.2, 0.3 (from bottom 
to top). The horizontal lines are tie-lines illustrating isotopic- 
nematic coexistence. 



Including quadratic terms leads to PcD^ = 1.243 in the 
C ^ limit. Because of the restricted number of allowed 
orientations, the I-N transition takes place at relatively low 
density compared to the predictions of simuLations for thin 
circular platelets with continuous orientationQiD. This point 
has already been discussed in the case of rod-like particles 
((^ > l)tj. The estimate ( |3^ ) predicts that pc increases upon 
increasing the aspect ratio C (< 1), in agreement-with recent 
experimental dataS and theoretical calculationsB. 

For C > 0, quantitatively more reliable phase diagrams 
may be derived by minimizing the grand potential (^) , with 
the FMT excess free energy density (|30|), with respect to 
the order parameter s for each density p. The resulting 
equations of state in the isotropic (s — 0) and nematic (s 7^ 
0) phases are plotted versus p in Fig. ||, together with the 
horizontal tie-lines connecting the two phases, for several 
aspect ratios C- The value of the effective packing fraction 
Vi = PiD^ and = pnD^ of the coexisting phase are 
seen to be rather insensitive to the aspect ratio Cj while 
the first order transition narrows as C increases, i.e., A77 = 
Vn — Vi decreases with increasing The reduced pressure 
P* — PD^ / {ksT), however, increases rapidly with as one 
might expect. 

In order to study the possible onset of an uniaxial nematic 
phase {px = Py), we rewrite the Euler-Lagrange equations 
in terms of p and s: 



In 



1 + 2s 2D{L - Dfps 



1-s 



1 - D^Lp 



Dfp's 



3{l-D^Lpy 



{D + 2L + [L - D]s) .(39) 



re- 



If only terms to cubic order in p are retained, Eq. ( |3£^ 
duces to Eq. (^7|). The bifurcation density follows from a 
low-s expansion of Eq. (B9h : 



G 



1 + GD^L 



(40) 



with 



G 




D^\{D + 2L)^ {D + 2L){L-Dy 



(41) 



Thus, for C ^ L/D ^ 0.2,0.4,0.6,0.8, p^''\ as cal- 
culated from Eq. (^l|), takes the values p^'^^D^ = 
1.259, 1.229, 1.165, 1.085, which are rather insensitive to the 
aspect ratio C- 



V. PHASE DIAGRAMS OF BINARY PLATELET 
MIXTURES 

We next consider the richer case of highly asymmetric 
binary mixtures of large and small platelets. Polydisper- 
sity, both in diameters and in thickness, is intrinsic to 
the expfidmental gibbsite samples which have been recently 
studiedBH. In the case of mixtures of spherical colloidal par- 
ticles large-asymmetry may lead to depletion-induced phase 
separationo. In the case of platelets depletion-induced seg- 
regation competes with nematic ordering, as is also the case 
for mixtures of long and short rodg£j, or of thin and thick 
rodsE3. Recent preliminary work on mixtures of continu- 
ously rotating thin circular platelets (C = 0) suggests that 
depletion-induced fractionation is-strongest when the large 
platelets are nematically orderecll3. Fraction|ation was also 
predicted for thin and thick circular plateletsQ, on the,basis 
of Parsons' heuristic extension of Onsager's theoryoES. 

In this section we consider binary mixtures of square 
platelets, using the FMT excess free energy (^0|). Various 
phase diagrams constructed as function of the chemical po- 
tential /i2 of the small platelets and the number density pi 
of the large platelets (Di > D2) are shown in Figs. |j, ||. 
The tie-lines are horizontal because of equality of p2 of the 
coexisting phases. For large negative chemical potential p2^ 
the systems exhibit a first oder I-N transition and the den- 
sity gap at the I-N transition broadens with increasing p2 
in agreement with our previous-calculations for binary mix- 
tures of infinitely thin plateletaHj. Upon increasing p2, the 
phase diagrams presented in Figs. ^, ^ differ qualitatively, 
and it is worthwhile to distinguish the following cases. 

Fig. 1^ (a), (b): The width of the I-N transition broadens 
continuously with increasing p2- At high densities in the 
nematic phase, the small platelets are nearly excluded and 
hence the large platelets are in equilibrium with a reservoir 
of small platelets. This implies that the pressure of the 
nematic phase is equal to that of a reservoir of small platelets 
at a chemical potential p2- 
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FIG. 4: Phase diagram of a fluid consisting of a binary mix- 
ture of small platelets [L2/D2 ~ 0.1] and large platelets [(a): 
D1/D2 = 10, Li/Di = 0.04; (b): D1/D2 = 5, Li/Di = 0.04; 
(c): D1/D2 ~ 10, Li/Di = 0.01] as a function of the chem- 
ical potential of the small platelets and the density of the 
large platelets pi. The straight solid lines are tie-lines illustrat- 
ing isotropic-nematic (I-N) coexistence. The dashed line in (c) 
marks coexistence of an isotropic phase (I) with two nematic 
phases (Ni and N2) of different compositions. 



Fig. ^ (c): Decreasing the aspect ratio (^1 = Li/Di 
(as compared to Fig. ^ (a), (b)) at fixed L2, D2 leads to 
a smaller width of the I-N transition and a shift of the 
I-N transition density to smaller values of pi as already 
discussed in the case of monodisperse platelet fluids (see 
Eq. (H)). The I-N transition ends in an I-N1-N2 triple 
point, at which two nematic phases (Ni, N2) coexist with 
an isotropic phase. Above the triple point there is coex- 
istence between a low-density isotropic and a high-density 
nematic phase. The nematic- nematic (N1-N2) coexistence 
region is bounded by a lower critical point. 

Fig. H (a), (b): The I-N coexistence regime is bounded 
by an upper critical point above which a single stable ne- 
matic phase is found. The nematic phase demixes into two 
nematic phases (Ni, N2) at sufflciently large values of p2- 

Fig. H (c): The width of the TN transition broadens with 
increasing aspect ratio ^2 = L2/D2 (as compared to Fig. ^ 
(a), (b)) at fixed Li, Di. Moreover, the lower critical point 
of the N1-N2 coexistence region shifts to smaller values of 
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FIG. 5: Phase diagram of a fluid consisting of a binary mixture 
of large thin platelets [L\/D\ — 0.01] and small platelets [(a): 
D1/D2 = 2, L2/D2 = 0.16; (b): D1/D2 = 1.66, L2/D2 = 0.13; 
(c): D1/D2 ~ 2, L2/D2 ~ 0.2] as a function of the chemi- 
cal potential of the small platelets fi2 and the density of the 
large platelets pi. The straight solid lines are tie-lines illustrat- 
ing isotropic-nematic (I-N) and nematic-nematic (N1-N2) coex- 
istence, respectively. The dashed line in (c) marks coexistence 
of an isotropic phase (I) with two nematic phases (Ni and N2). 



H2, until the N1-N2 and I-N coexistence regions start to 
overlap, giving rise to a triple point. 

A simple estimate of the critical aspect ratios required 
for a depletion-driven demixing transition is obtained from 
a second virial approximation of the excess free energy func- 
tional (16). A straightforward analysis shows that demixing 
at constant volume is possible, provided the following inter- 
action parameter x(si, S2) > 0: 



Xisi,S2) = '2Ei2{si,S2)-Eii{si,Si)-E22{s2,S2), (42) 

with 

Eij , 5j ) 

i(A + D^) [2{D,Dj + L,;L,)(1 - 

+ (Aij + DjLi){l + 2s,sj) + i{D,L, + DjLj)] ,(43) 

where i,j — 1,2. The values of the nematic order param- 
eters si and S2 are obtained by a numerical minimization 




0.00 
0.00 



FIG. 6: Phase diagram of a binary platelet fluid, corresponding 
to Fig. ^ (a), in the density-density, i.e., pi-p2 plane. The dotted 
lines indicate phase boundaries. The chemical potential of the 
small platelets is kept fixed for each solid line and increases from 
bottom to top: ^2/(fcsr) = -5, -4, -3, -2, -1, 0, 1, 1.9, 3. 



of the grand potential. In order study a possible nematic- 
nematic demixing transition one may investigate the inter- 
action parameter x(l, 1): 



x(i,i) 

Dl 



D2 



C1+C2 



D2 



8C1 - 8C2 



D2 
Dl 



, (44) 



where Ci — Li/Di and C2 = L2/D2- For the model param- 
eters used in Fig. ||, x(l, 1) > 0, in agreement with that fact 
that nematic-nematic demixing has been found by numeri- 
cal minimization of the FMT grand potential. 

Figure 6 displays an alternative representation of the 
phase diagram shown in Fig. ^ (a) in terms of the num- 
ber densities of the large and the small platelets pi and p2, 
respectively. The figure illustrates how the composition of 
the mixture varies upon increasing the chemical potential 
(and hence the reservoir density) of the small platelets. 



VI. CONCLUSIONS 

We have generalized Zwanzig's model of cuboids (or par- 
allelepipeds) with only three allowed orientations to the 
case of monodispersc and bidisperse systems of platelets 
(C = L/D < 1). In view of the uncertainties concerning 
the convergence of the virial series, we have used Rosen- 
feld's FMT|V-|along the lines of Cuesta's work on parallel 
hard cubesHH to derive a free energy function of the den- 
sity and the nematic order parameter, which on the basis 
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of earlier experience with other hard core systems, is ex- 
pected to be quantitatively accurate. In the monodisperse 
case, the isotropic to nematic transition narrows with in- 
creasing aspect ratio C. Very rich phase diagrams, involv- 
ing an isotropic and one or two nematic phases of different 
concentrations are found for bidisperse systems; the phase 
diagrams are very sensitive to the size and aspect ratios and 
although our results can only be considered to be of quali- 
tative significance, due to the restricted number of allowed 
orientations, we believe that these results clearly point to 
the possibility of tuning the phase behaviour by appropri- 
ate choices of the relative diameters and thicknesses of real 
platelets. This is of genuine technological importance for 
the design of clay based drilling fluidil. Since clay platelets 
are highly charged, Coulombic interactions must be included 
in the model and this could lead to major changes in the 
predicted phase diagrams. Calculations of binary platelet 
phase diagrams, including Coulombic forces, along the lines 
of Ref.Ej are under way. 
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APPENDIX A: THIRD AND FOURTH VIRIAL 
COEFFICIENTS FOR THE ZWANZIG MODEL OF 
PLATELETS 

According to the standard definition of the third virial co- 
efRcient involving three particles of species (3i, (32, (where 
here Pi = x,y ot z denote one of of the three possible ori- 
entations of the platelets): 

i?3(/3l,/32,/53) = 

drdr' //3i,/32(r)//3i,/33(r')//32,/33 (i" " r') 

00 00 

ijl J dra J dr'^e{Sa;i3^,i32 -\ra\) 

^ —00 —00 

xBiSc^.p^f,, ~ \r'J)eiSa;M, ~ |r„ ~ r'J) . (Al) 



where Sa-^p-i^^jS^ = {Sa;i3i + Sa-iSz)/^ and definitions (||) and 
(^) have been used. The two-dimensional integrals involving 
Heaviside step functions are easily calculated, and lead to 
the following results: 



-83(2;, X, x) 
B3{x,x,y) 



dL^D"^ , (A2) 
LD^{L + 2D){D + 2L), (A3) 

B3(x,y,z) = ^D^2L + D)\ (A4) 

while for the isotropic phase (where px = Py = Pz ^ p/'^), 
B3, as given by (|9|), reduces to: 



16 



47 



52 



B3 = —D" + —D^L + —D^L^ + —D^L 



3t3 



27 



27 



(A5) 



These expressions are easily generalized to the case of binary 
mixtures of platelets, with the added complication of species 
indices. 

Similarly, there are 4 independent 4-th virial coeffi- 
cients in the one-component case, corresponding to vari- 
ous orientations of the four particles, namely B4{x, x, x, x), 
B4{x,x,x,y), Bi{x,x,y,y) and B4^{x,x,y, z). Each of 
these has conlxibutions from 4 irreducible diagrams (see 
Eq. (4.5.13) iro) for each set of orientations. The calcu- 
lations of the multiple integrals are tedious, and we have 
evaluated B^ only for three limiting cases. We find for the 
isotropic phase: 

(a) Long rods, corresponding to the limit C — > 00: 



Ba 



32 
243' 



3t& 



-D-'L' 



(A6) 



(b) Thin platelets, corresponding to the limit C ^ 

B, = -^D^. 
243 

(c) Parallel cubes, corresponding to C = 1: 

34 



(A7) 



Bi 



(A8) 



In the case of long rods and thin platelets the fully con- 
nected diagram does not contribute to B4, and the calcula- 
tions of the other diagrams is considerably simplified. 
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